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We propose a new mechanism for turbulent mean-field dynamo in which the magnetic fiuctuations resulting 
from a small-scale dynamo drive the generation of large-scale magnetic fields. This is in stark contrast to the 
common idea that small-scale magnetic fields should be harmful to large-scale dynamo action. These dynamos 
occur in the presence of large-scale velocity shear and do not require net helicity, resulting from off-diagonal 
components of the turbulent resistivity tensor as the magnetic analogue of the “shear-current” effect. Given the 
inevitable existence of non-helical small-scale magnetic fields in turbulent plasmas, as well as the generic nature 
of velocity shear, the suggested mechanism may help explain generation of large-scale magnetic fields across a 
wide range of astrophysical objects. 


Astrophysical magnetic fields are observed to be well- 
correlated over length and time scales far exceeding that of 
the underlying fluid motions. Beautiful in its regularity, the 
22 -year solar cycle is the most well-known example of this 
behavior m. Such large-scale structure is puzzling given that 
strong magnetic fields are expected to emerge through the 
stretching and twisting of field lines by smaller scale turbu¬ 
lence. As the primary theoretical framework to study such be¬ 
havior, mean-field dynamo theory examines how large-scale 
magnetic fields develop due to these small-scale turbulent mo¬ 
tions. This splitting between scales is captured by the mean- 
field average; the average of a fluctuating quantity vanishes by 
definition (b = 0 ), while the average of a large-scale field is 
itself (B = B). An average of the induction equation, which 
governs evolution of the magnetic field within magnetohydro¬ 
dynamics (MHD), leads to lO 

d,B = Vx{UxB) + VxS+—V'^B, (1) 

Rm 

where Rm is the magnetic Reynolds number, a dimensionless 
measure of the plasma resistivity, and U and B are the large- 
scale velocity and magnetic field. The electromotive force, 
S = u X b,is 3.n average of the fluctuating fields (u and b) and 
responsible for dynamo action. In the early phases of a dy¬ 
namo, the mean fields can be considered a small perturbation 
to the underlying turbulence. Combined with an assumption 
of scale-separation between small-scale and mean fields, this 
allows a Taylor expansion |[3l[4l of 6 in terms of B, 

S = aoB-\-/3o VB + • • • , (2) 

where a, fS are the tensorial transport coefficients, calculated 
from the small-scale fields m Since these depend on the 
large-scale fields, a solution to Eq. O requires knowledge of 
how S changes with B (and possibly U), essentially a statisti¬ 
cal closure for inhomogenous MHD. 

Historically, much work has focused on kinematic dynamo 
theory, in which u is uninfluenced by the magnetic field 
Eia. Kinematic theory predicts large-scale dynamo instabil¬ 
ity when the fluid motions possess helicity, / m • V x m Jjc 9 ^ 0. 
However, the applicability of such predictions has been called 
into question by a number of authors In particular, 


above modest Reynolds numbers in both helical and non¬ 
helical flows, the small-scale dynamo HI causes b to grow and 
saturate much more rapidly m than B. This violates the kine¬ 
matic assumption, both because u is altered before B grows 
significantly, and because a dynamically important b exists 
independently of B. The buildup of small-scale fields is the 
origin of “a quenching”, in which the mean field saturates 
well before reaching amplitudes consistent with observation 
ifTomsl due to the adverse influence of b. 

In this letter we show that in turbulence with large-scale 
velocity shear, it is possible and realizable to have the small- 
scale dynamo enhance the growth of the large-scale dynamo. 
We demonstrate this both with statistical simulation ina, in 
which the effect is very clear but applies rigorously only at 
low Reynolds numbers, and through calculation of transport 
coefficients from direct numerical simulations (DNS). In ad¬ 
dition, the existence of the effect has been confirmed analyt¬ 
ically using the second-order correlation approximation El, 
which agrees with previous spectral r approximation calcula¬ 
tions El. 

All computations presented here use the incompressible 
MHD model in the shearing box, employing homogenous 
Cartesian geometry and periodic boundary conditions in the 
shearing frame. With a mean flow Uq = -S xy imposed across 
the domain, this setup is designed to represent a small “patch” 
of turbulent fluid in large-scale velocity shear. We force with 
non-helical white-in-time noise at small scales and study the 
generation of larger-scale magnetic fields, in a way similar to 
previous authors EKTH. The mean-field average is defined 
as an average over the horizontal (x and y) directions, such 
that the mean magnetic fields B depend only on z. We also 
allow for system rotation through a mean Coriolis force, since 
shear typically arises due to differential rotation in astrophys¬ 
ical objects. The rotation O is aligned with z (antiparallel to 
V X t/o), perpendicular to the flow Uq. 

For the chosen horizontal average, inserting Eq. ^ into Q 
gives 

~ ~^yX^Z^X ~ Otyyd^By — TJyxd^By -f Tjfyd^Bx 

dtBy = -SBx + axxdzBx + otxyd^By - Pxy^lBx + Oix^lBy, (3) 
using velocity shear Uq but neglecting other mean velocities, 
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and defining tjh = riu + Rm“^ Here atj and rjij are the a 
effect and turbulent resistivity tensors respectively, with the 
4 components of T]ij relatable to the elements of the full 
tensor [Eq. Due to homogeneity and reflectional sym¬ 
metry (vanishing net helicity), atj must vanish when averaged 
over a suitably large time or number of realizations m, and 
indeed our measurements confirm this. There is no such con¬ 
straint on rjij, and rjyx is very important throughout this work 
due to its coupling with the shear. In particular, neglecting 
fluctuations in a and assuming diagonal resistivities are equal 
= rjty = rft), the least stable eigenmode of Eq. ^ for a 
mode of vertical wavenumber k grows at 

y = k [T]yx(-S + k^T]xy)f^^ - k^TJt. ( 4 ) 

Since S » rfij, dynamo action is possible without an a effect 
if riyx < 0. 

Subsequent to early analytic work HMD, it was found 
kinematically that rjyx > 0 (at least at low Rm), and several au¬ 
thors have thus concluded that a coherent shear dynamo can¬ 
not explain observed field generation |[T3[22l|23l. Instead, a 
popular theory is that temporal fluctuations in atj cause an in¬ 
coherent mean-field dynamo. Importantly, in such a dynamo, 
B (z, t) cannot have a constant phase in time as it grows, since 
the average of B over an ensemble of realizations vanishes, 
implying B must be uncorrelated with itself after t > 

1241 . While incoherent dynamos are possible in a variety of 
situations, here we argue for a different situation - magnetic 
fluctuations act to substantially decrease and potentially re¬ 
verse the sign of riyx, causing the onset of a coherent large- 
scale dynamo that can overwhelm the incoherent dynamo. 

Our first method illustrating this effect is quasi-linear statis¬ 
tical simulation; see Ref. Il26l for further details. The method 
starts by forming equations for the fluctuating fields u and b, 
and linearizing these; i.e., neglecting fluctuation-fluctuation 
nonlinearities such as m • Vm and u Vb. One can then derive 
an equation for the fluctuation statistics, C = iXiXj) (where 
X = (u,b) and (•) denotes the average over an ensemble of re¬ 
alizations), as a function of U and B |[T4l[27l[28l. Einally, us¬ 
ing 8x = C 26 - C53 and 8y = C34 -Cei, 8 can be fed directly 
into Eq. 0, resulting in a closed system of equations. Note 
that the method does not assume a scale-separation between 
mean and fluctuating fields. Importantly, since the statistics 
are calculated directly, an incoherent dynamo is not possible, 
and statistical simulation offers a direct probe of the coherent 
effect. The linearity of the fluctuation equations eliminates 
the small-scale dynamo; accordingly, to excite homogenous 
kinetic and magnetic fluctuations, both u and b are forced 
at small scales (kf = 671) with the statistics of white-in-time 
noise. The resulting MHD turbulent bath could be considered 
as some approximation to kinetically forced turbulence after 
saturation of the small-scale dynamo. 

To study the magnetically driven dynamo, we keep the total 
forcing level constant, successively increasing the proportion 
of magnetic forcing from purely kinetic 1^ . Results with 
shear but no rotation are illustrated in Eig.[2 The presence of 



FIG. 1. Time development of the mean-field energy, Eb - 
f dzBY2, in quasi-linear statistical simulation. Small-scale fields 
are forced at kf - Gn, Rm = Urmsirikf ^ 5 (here 77 is the resis¬ 
tivity, Pm = Rm/Re = 1), S' -2 and the box has dimensions 
[Lx^Ly,L^ - (1,1,4) with resolution (28,28,128). As well as forc¬ 
ing the momentum equation, the induction equation is forced to ex¬ 
cite homogenous magnetic fluctuations, emulating a small-scale dy¬ 
namo. Total forcing is kept constant in each simulation but the pro¬ 
portion of magnetic forcing (fblf) is increased from 0 to 0.8 (25). 
As fb! f increases, the growth rate of the mean field increases due to 
the change in sign and subsequent decrease of riy^. 


the magnetically driven dynamo is evident, becoming slightly 
unstable when magnetic forcing accounts for 0.4 of the total 
and increasing the growth rate thereafter. This sustained pe¬ 
riod of exponential growth due to magnetic fluctuations is not 
possible to see in DNS, since the mean field will immediately 
come into approximate equipartition with the small-scale field 
due to the finite size of the system. 

The formal applicability of statistical simulation is limited 
to low Rm due to the quasi-linear approximation. Our sec¬ 
ond method utilizes DNS, forced kinetically at small scales, 
to show that magnetic fields arising consistently through the 
small-scale dynamo can drive a coherent large-scale dynamo. 
To this end, we directly calculate transport coefficients from 
nonlinear simulation before and after the saturation of the 
small-scale dynamo. We use the incompressible MHD code 
SNOOPY 1291 . which uses the pseudo-spectral method and 
shearing periodic boundary conditions. The chosen Reynolds 
numbers fm are moderate - small enough such that there 
is no self-sustaining turbulence in the absence of small-scale 
forcing (although effects may be similar even when this oc¬ 
curs (sol) - and ensembles of 100 simulations are run with 
shear Uq = Sxy, both with and without Keplerian rotation. 
At these parameters, the prevalence of the coherent large-scale 
dynamo depends on the realization (see Eig.[^, and it appears 
that the coherent effect cannot always overcome fluctuations 
in 8 immediately after small-scale saturation, although the dy¬ 
namo develops after a sufficiently long time [e.g.. Pig. |^d) 
near t = 150]. This behavior seems generic when the coherent 
dynamo is close to its threshold for excitation and we have ob¬ 
served similar structures when the induction equation is driven 
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FIG. 2. Example spatiotemporal By evolutions for non-rotating (a-b) 
and Keplerian rotating (c-d) turbulence at Rm = Urmsh^f ~ 15 {kf = 
6n, Tj = 1/2000, Pm = 8), 5 = 1, in a box of dimension (1,4,2) with 
resolution (64,128,128). The first examples in each case [(a) and 
(c)] show By when a coherent dynamo develops, while the second 
examples [(b) and (d)] illustrate the case when it is more incoherent. 
The main factors in distinguishing these are the coherency in phase 
of By over some time period and the amplitude at saturation, which is 
larger in the coherent cases. In general, the rotating simulations are 
substantially more coherent. The hatched area illustrates the region 
of small-scale dynamo growth. The fitting method used to compute 
transport coefficients (see Fig.[^ is applied between the dashed lines 
{t = 50^ 100). 
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FIG. 3. Measurements of the turbulent transport coefficients for 
100 realizations of the simulations at the same parameters as those 
in Fig.[^ (a) rjxx and a^y coefficients (x and y axes respectively), no 
rotation, (b) T]yx and ayy coefficients, no rotation, (c) j]xx and a^y co¬ 
efficients, rotating, (b) T]yx and ayy coefficients, rotating [see Eq. |^]. 
Unfilled markers in each plot show coefficients measured from each 
of the individual realizations, with mean values displayed with solid 
markers and the shaded regions indicating error in the mean (2 stan¬ 
dard deviations). Black markers illustrate the kinematic transport 
coefficients, with grey shaded regions indicating the error. After sat¬ 
uration of the small-scale dynamo, we calculate aij(t) and T]ij(t) by 
solving Eq. ^ approximately (see text), taking the mean from t = 50 
to t - 100. This limited time window is chosen to avoid capturing 
the saturation phase of the large-scale dynamo, since rjij is presum¬ 
ably modified in this phase. In both methods of computing transport 
coefficients, a coefficients are zero to within error as expected, and 
the scatter between simulations is of a similar magnitude to that of 
rjij if one accounts for their different units (it is necessary to divide a 
by a characteristic k value). 


directly at lower Rm 1261 . Notwithstanding this variability in 
the dynamo’s qualitative behavior, measurement of the trans¬ 
port coefficients illustrates that the rfy^ coefficient decreases 
after the magnetic fluctuations reach approximate equiparti- 
tion with velocity fluctuations at small scales. 

At low times, we use the test-field method to measure the 
kinematic a and i], fixing the mean field and calculating S, 
with no Lorentz force ^ [TSl. Since the small-scale dy¬ 
namo grows quickly, test-fields are reset every t = 5. After 
small-scale saturation, standard test-field methods are inap¬ 
plicable ll3n . Instead, we extract B and S simulation data 
and calculate Tjij(t)) directly from Eq. ^ by comput¬ 

ing f dzSiS for each of 0 = (Bx^By^d^Bx^d^By) and solv¬ 
ing the resulting matrix equations (in the least-squares sense) 
at each time point. This method is similar to that presented 
in Ref. ; however, we additionally impose the constraints 
ByyiO ~ dxxif)') ^xxif) — ^yyif) 3.nd (Xyxit) = dxyif) — A- While 
these changes may appear to make the method less accurate, 
they in fact achieve the opposite by reducing the influence of 
Bx. This is necessary because Bx has a high level of noise in 


comparison to its magnitude, and because this noise is cor¬ 
related with the noise in By (due to Bx directly driving By) 
and &y (through dtBx = -dzSy -f • • • )• These correlations are 
very harmful to the quality of the fit, causing unphysical neg¬ 
ative values for rjyy 1 ^ . which also pollute measurement of 
other coefficients. It is straightforward to show that the sys¬ 
tematic errors caused by our constraints on the transport coef¬ 
ficients are less than 1% for the dynamos in Fig.[^ so long as 
dxx ~ dyy when time averaged. We have verified the method is 
accurate by comparison with the kinematic test-field method 
at low Rm, where there is no small-scale dynamo |[2^[33l and 
the rotation can be used to explore a range of r]yx. In addi¬ 
tion, the measurements are independently verified in Fig. 
below oa Due to the short time-window, measurements 
of the transport coefficients after small-scale saturation vary 
significantly between realization, as should be expected from 
Fig.[^ Nonetheless, an average over the ensemble of 100 sim¬ 
ulations illustrates a statistically significant change in rjyx that 
is consistent with observed behavior. 
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FIG. 4. Evolution of the mean-field magnitude for a sample of the 
ensemble of rotating simulations (Figs.j2jj^. Here 5, the mean-field 
magnitude, is i\B\^ + where ^is the largest scale Fourier 

mode of Bf. In each plot the solid blue curve shows data taken from 
the simulation. The dashed red curve shows the corresponding ex¬ 
pected evolution, using the time-dependent calculated values of the 
transport coefficients, smoothed in time using a Gaussian filter of 
width 5. Finally, the dotted black curve illustrates the expected evo¬ 
lution with all a coefficients artificially set to zero. The similarity 
between this evolution and that including a (dashed curve) illustrates 
that in many (but not all) cases the dynamo is primary driven by rfy^. 
We list the measured mean of in each plot to show that lower val¬ 
ues generally lead to stronger mean-field growth, as expected for a 
coherent dynamo. For reference, at the measured rj^x ~ 0.006, the 
coherent dynamo is unstable below rjyx - -0.00036. 


We thus conclude that small-scale magnetic fluctuations act 
to decrease rjyx, and that in some realizations (or after a suf¬ 
ficiently long time period) a coherent large-scale dynamo de¬ 
velops as a result. 

To summarize, in this letter we have demonstrated that 
small-scale magnetic fluctuations, excited by small-scale dy¬ 
namo action, can drive large-scale magnetic field generation. 
The mechanism is a magnetic analogue of the “shear-current” 
effect arising through the off-diagonal turbulent re¬ 

sistivity in the presence of large-scale shear flow. We have 
demonstrated its existence numerically using both DNS, with 
measurements of mean-field transport coefficients before and 
after small-scale dynamo saturation, and through quasi-linear 
statistical simulation. 

More work is needed to precisely assess regimes in which 
the magnetically driven dynamo might dominate, as well as 
its behavior at higher Reynolds numbers where self-sustained 
turbulence is possible Ool. Another interesting question re¬ 
gards whether a magnetic dynamo can remain influential in 
the presence of net helicity and an a effect, particularly as 
small-scale dynamo may be suppressed by shear fSSl . While 
such questions may be difficult to answer definitively, the 
generic presence of magnetic fluctuations in plasma turbu¬ 
lence gives us some confidence that the proposed mechanism 
could cause large-scale dynamo growth in the wide variety of 
astrophysical systems with velocity shear. 

This work was supported by a Procter Fellowship at Prince¬ 
ton University, and the US Department of Energy Grant 
DE-AC02-09-CH11466. The authors would like to thank 
A. Schekochihin, J. Krommes, and 1. Rogachevskii for en¬ 
lightening discussion and useful suggestions. 


Results are illustrated in Eig. In the kinematic phase 
without rotation, we see rjyx = (4.1 + 1.6) x 10“"^ in qual¬ 
itative agreement with previous studies CD. With rotation 
rjyx = (0.6 + 1.2) X 10“^, consistent with a reduction in rjy^ 
due to the n X / effect 0. After saturation of the small-scale 
dynamo, rjyx = (-0.1 + 1.0) x 10“"^ for the non-rotating case, 
while T]yx ~ - (2.0 + 0.8) x 10“"^ in the rotating case - the 
same reduction in each to within error. Values for the diag¬ 
onal resistivity are smaller after saturation, as expected since 
the velocity fluctuation energy decreases (by a factor ~ 1.4). 
The values of (rjxx, rjyx) show that the dynamo is slightly stable 
on average in the non-rotating case and marginal in the rotat¬ 
ing case. However, the coefficients vary significantly between 
realizations, sometimes yielding larger growth rates, and mea¬ 
surements match observed growth of the mean field for indi¬ 
vidual realizations. We illustrate this in Eig.|^ which demon¬ 
strates consistency between the measured transport coeffi¬ 
cients and observed mean-field evolution by solving Eq. 0 
directly [using the time-dependent Tjij(t) and for a sam¬ 

ple of the rotating simulations. In addition, by artificially re¬ 
moving we illustrate that cases with more negative rjyx 

are driven primarily by this, rather than a stochastic-a effect. 
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